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This  report  is  the  final  report  on  Contract  No.  F44620-74-C-0057  to 
study  the  application  of  ultrasonic  elastic  wave  scattering  to  the  detection, 
identification  and  characterization  of  defects  in  solids.  The  contract  was 
initiated  on  February  1,  1974,  and  terminated  on  March  15,  1977.  During  the 
three-year  period  it  entailed  approximately  1-1/2  man-years  of  effort  of 
which  approximately  one-third  was  experimental  measurement  and  two-thirds 
theoretical  analysis.  During  the  first  two  years  the  theoretical  analysis 
was  compared  with  experimental  measurements,  but  it  became  apparent  that  this 
approach  was  not  adequately  providing  a validation  of  the  calculations  since 
any  discrepancies  could  be  ascribed  to  various  experimental  shortcomings  such 
as  instrumental  sensitivity  and  angular  resolution,  lack  of  monochromaticity 
in  the  transducer  pulses,  variability  in  the  physical  parameters  of  the 
titanium  alloys,  and  attenuation  and  scattering  of  the  ultrasonic  signals  by 
metallurgical  inhomogeneities  (second  phases) . For  this  reason  the  effort  of 
the  third  year  of  the  program  was  entirely  analytic. 
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II.  SUMMARY  OF  PROGRESS 

The  results  of  the  first  two  years  have  been  reported  in  two  annual 
reports  [Cohen,  1975;  Cohen,  1976]  and  an  interim  report  [Cohen  and  Tittmann, 
1975]  and  will  only  be  summarized  here.  For  details,  the  reader  is  referred 
to  the  prior  reports  which  must  be  considered  as  part  of  the  full  documenta- 
tion of  this  contract.  The  theoretical  results  on  the  formulation  of  the 
scattering  of  elastic  waves  will  be  reported  here  in  some  detail. 

A.  Scattering  of  Elastic  Waves  by  Spherical  Cavities  and  Inclusions 

In  order  to  verify  the  experimental  techniques  and  to  develop  the  necessary 
theoretical  basis,  the  initial  task  of  this  contract  was  the  review  of  the 
existing  calculations  of  scattering  by  spheres  [Cohen,  1975]. 

B.  Scattering  of  Elastic  Waves  by  Spherical  Cavities 

Spceimens  of  titanium  alloy  containing  a spherical  void  of  diameter 
0.4  mm  or  0.8  mm  were  prepared  by  diffusion  bonding  techniques.  Experimental 
data  on  the  scattering  of  ultrasonic  pulses  by  these  cavities  were  compared 
with  the  theoretical  calculations.  The  data  were  in  reasonable  agreement  with 
theory.  Discrepancies  could  be  ascribed  to  the  lack  of  monochromaticity  of 
the  pulses  at  high  frequencies,  and  to  the  scattering  and  attenuation  of  the 
pulses  by  the  metallurgical  inhomogeneities  in  the  titanium  alloy  [Cohen,  1975]. 

C.  Scattering  of  Compressional  and  Shear  Waves  by  a Spherical  Inclusion 

A tungsten  carbide  inclusion  (commercially-available  tungsten  carbide 
ball,  nominal  diameter,  1/32",  0.794  mm)  was  inserted  into  a diffusion-bonded 
titanium  alloy  block.  Measurements  of  compressional  wave  scattering  were 
carried  out  at  nominal  frequencies  of  2.25,  5.0  and  10.0  MHz.  The  transducers 
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used  had  low  Q (approximately  5)  and  hence  not  adequately  monochromatic. 

Because  of  this,  waveform  synthesis  was  required.  The  2.25  MHz  transducer 
showed  a characteristic  with  a full  width  at  half  maximum  of  0.8  MHz.  The 
5 MHz  transducers  had  a peak  output  at  4.5  MHz  with  a full  width  at  half 
maximum  of  1.9  MHz,  and  the  10  MHz  transducers  peaked  at  7.4  MHz  with  a full 
width  at  half  maximum  of  2.6  MHz. 

Experiments  with  incident  shear  waves  were  carried  out  with  high  Q shear- 
cut  quartz  plates  bonded  to  the  faces  of  the  sample. 

These  experiments  were  in  agreement  with  theoretical  calculations  of  the 
cross  sections  except  for  the  amplitude  of  cross-polarized  (depolarized)  shear 
wave  scattering  at  5 MHz.  The  precise  reasons  for  the  anomalous  appearance  of 
this  depolarized  amplitude  (an  obversed  scattered  shear  wave  polarized  at 
right  angles  to  the  scattering  plane  when  the  incident  shear  wave  is  parallel 
to  that  plane,  which  is  forbidden  by  symmetry  for  scattering  by  a sphere)  is 
not  understood  but  has  been  ascribed  to  the  bulk  scattering  by  the  metallurgical 
second  phase  in  the  titanium  alloy  [Cohen,  1976]. 

D.  Reformulation  of  the  Calculation  of  Scattering  by  Spheres 

The  analyses  of  Ying  and  Truel  [Ying  and  Truel,  1956]  on  the  scattering  of 
compressional  waves  by  spherical  scatterers  and  of  Einspruch,  Witterholt  and 
Truel  [Einspruch,  et  al. , 1960]  on  the  scattering  of  shear  waves,  were  rederived 
in  a unified  format  which  makes  possible  several  useful  extensions  [Cohen  and 
Tittmann,  1976]: 

a.  Both  scattering  processes  are  expressed  in  terms  of  the  same  basic 
"scattering  matrix"  which  allows  the  amplitudes  of  the  spherical 
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harmonics  of  the  scattered  wave  field  to  be  expressed  in  terms  of 
those  of  the  incident  wave  field. 

b.  This  same  matrix  applies  also  to  the  scattering  of  any  incident  wave 
field  (not  necessarily  plane)  and  hence  is  in  a convenient  form  to 
permit  the  analysis  of  near-field  scattering  or  the  multiple 
scattering  of  two  spheres. 

c.  The  reciprocity  theorem  for  mode  conversion  scattering  by  a sphere: 

At^<9)  ■ kVt«> 

may  be  established  directly  from  the  symmetry  of  the  "scattering 
matrix. " 

E.  "Corrected"  Born  Approximation 

The  Born  approximation  is  based  on  an  integral  equation  formulation  of 
the  scattering  process  (see  Section  III.D  below).  The  total  displacement  field 
is  expressed  as  a sum  of  the  imposed  incident  field  plus  an  integral  over  the 
volume  of  the  scatterer  which  involves  the  displacement  field  in  the  scatterer 
and  the  differences  between  the  material  constants  of  the  scatterer  and  the 
constants  of  the  external  host  medium.  The  first  Born  approximation  replaces 
the  displacement  field  in  this  integral  by  the  incident  displacement  field. 

The  resultant  expression  for  the  scattering  cross  section  is  an  adequate 
approximation  only  when  the  properties  of  the  scatterer  differ  slightly  from 
those  of  the  host  medium.  A correction  to  this  first  order  Born  approximation 
has  been  developed  [Cohen,  1976]  which  yields  the  exact  long-wave  length  limit 
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for  the  scattering  cross  sections  for  spherical  scatterers  for  arbitrary  varia- 
tion of  scatterer  parameters.  In  particular,  this  corrected  Born  approximation 
gives  the  correct  foward  scattering  amplitude  for  a spherical  cavity  in  the 
Rayleigh  limit. 

A further  extension  of  the  Born  approximation,  which  is  developed  in 
Section  III.E  below,  gives  results  which  are  exact  in  the  Rayleigh  limit  for 
scatterers  of  arbitrary  shape  (subject  at  most  to  the  restriction  of  a smooth 
boundary) . 

F.  Variational  Formulation 

From  a variational  formulation  of  the  cross  section,  a further  extension 
of  the  "corrected"  Born  approximation  is  possible.  This  includes  the  long- 
wavelength  extension  mentioned  in  the  previous  paragraph  for  cj  “*■  0.  The  range 
of  its  applicability  as  a function  of  u)  (or  of  ka,  <a)  may  be  partially 
evaluated  by  comparison  with  the  known  exact  solution  for  a sphere.  Applica- 
tion of  this  formulation  to  ellipsoids,  in  particular  to  ellipsoids  of  revolu- 
tion, is  underway.  Further  extensions  which  should  increase  the  range  of 
applicability  of  this  class  of  approximation  are  possible  without  significantly 
increasing  the  complexity  of  the  calculational  algorithms.  Details  are  given 
in  Section  III.F. 
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III.  DETAILS  OF  THIRD  YEAR'S  RESULTS— GREEN' S FUNCTION, 
RECIPROCITY  AND  THE  BORN  APPROXIMATION  FOR  THE 
SCATTERING  OF  ELASTIC  WAVES 

A.  Differential  Equation  and  Green's  Function 

The  displacement  field  in  an  elastic  medium  satisfies  the  equation 

pu  (r)  = (c  .u  . ) + f ( 

m — \ mspj  p,j /,s  m 

where  f is  the  applied  force  field.  In  Eq.  (1)  the  elasticity  tensor 

c . (r)  may  be  a function  of  position.  For  an  isotropic  medium  c . is 
mspj  — mspj 

given  by 


c , = \6  5 . + 

mpsj  ms  pj 


U(S  6 . 
tnp  sj 


<5  .6  ) 

mj  sp 


(2) 


and  for  any  arbitrary  medium  must  satisfy  at  least  the  symmetry  with 

respect  to  interchange  of  indices  expressed  in  Eq.  (2);  i.e., 

c . = c = c , = c , 

mspj  smpj  msjp  pjms 

For  displacements  which  are  harmonic  in  time  with  frequency  u)  Eq.  (1) 
becomes 


(c  .u  . ) 
v mspj  p,  j ' ,i 


2 

poo  u 


f = 
m 


(3) 


The  Green's  function  for  the  system  is  defined  as  the  displacement  field  pro- 
duced by  a delta  function  impressed  force  at  the  point  r , and  hence  GU(r;r  ) 

o p o 

is  the  solution  of  the  equation 
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(c  ,Gn  . (r;r  ))  + pto2Gn(r;r  ) + 6 6( r-r  ) 

' mspj  p,j o^s  m o mn o 


The  Green's  function  can  in  principle  be  defined  for  an  inhomogeneous  aniso- 
tropic system  as  the  solution  of  Eq.  (4),  but  it  is  useful  for  our  purpose  of 
exploring  the  scattering  of  elastic  waves  to  restrict  the  Green's  function  to 
the  solution  of  Eq.  (4)  for  an  infinite  homogeneous  isotropic  medium,  and 
hence  to  consider  only  the  Green's  function  which  is  the  solution  of 


(A+;j)Gn  (r;r  ) + uGn  (r;r  ) + ptj2Gn(r;r  ) + 5 5 (r-r  ) = 0 (5) 

o,pm o m,pp  — — o m o mn o 


with  the  boundary  condition 


lim  G (r;r  ) = 0 

I | _ m o 

r r* 


It  is  clear  that  the  solution  of  Eqs.  (5)  and  (6)  possess  translational 
invariance  and  reciprocity: 


<£< r;r  ) = G“(|r-r  |)  = G“(|r  -r | ) , 

m — — o m o n — o 


The  Green's  function  can  be  explicitly  written  [Morse  and  Feshbach,  1953; 
Gubernatis,  et  al. , 1976] 


G"(  I r-r  ) = G.  (R)  = 
m o jm 


1 

r 0 -i<R 

2 x e 

32 

-ikR 

e 

' 

, 2 
4tt  pcu 

< 0 . ^ i 

jm  R 

3x . 3x 
3 m 

( R 

R 

J 

We  now  introduce  the  spherical  Bessel  function  hQ(z)  and  write 
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T 2 

G.  (R)  = K36.  h (kR)  + (kcR)  - kh  (kR)  [ 

im  . . „ 2 jm  o ox . ox  I o o 1 

J 4Tripo)  [_  J 2 m ' 


where  < and  k are  the  wave  numbers  for  shear  and  compressional  waves. 


2 2 2 2 
k = pw  /y  , k = poj  /X+2\i 


hQ(z)  = ie  /z  = jQ(z)  - iyQ(z) 


We  shall  also  make  use  of  the  Bessel  function  relationships 


h , (z ) + h (z ) = — ■ h (z ) 


2 2 2 2 
and,  with  R = x^  + X£  + x^, 


a hn(kR)  . WkR) 
9xj  Rn  ' ^ Rn+1 


Using  these  expressions,  we  may  write  Eq.  (9)  in  the  form 


G.  (R)  = — w 6.  ^2K3h  (<R)  + k3h  (kR)z 

Jm  12it  ipu>2  0 0 > 


+ [~T-6 -k3h2*R>i- 
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The  symmetry  and  reciprocity  relationships  of  Eq . (7)  are  immediately 
evident  in  the  explicit  form  of  Eq.  (8)  or  Eq.  (14);  however,  a much  more 
general  reciprocity  expression  may  be  developed  directly  from  Eq.  (4)  for 
the  general  case  [Knopoff  and  Gangi , 1959;  Cohen,  1976].  We. rewrite  Eq.  (4) 
with  a slightly  different  notation 


(CmspjUp,j  (— 5— 1^), s + Pw2um(^l}  + = 0 


(15.1) 


V 

where  u^rjr^)  is  the  displacement  ar  _r  resulting  from  a unit  force  in  the 


v-direction  at  ;r^,  and 


(CmspjUp,j  2)),s  + + 5 fi(r-r2)  = 0 . (15.2) 


If  (15.1)  is  multiplied  by  u^(r;r  ) and  (15.2)  is  multiplied  by  uV(r;r  ) and 

m 2 m — — 1 

the  two  expressions  subtracted,  one  obtains 
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Now,  from  the  symmetry  of  the  elastic  coefficients  the  final  integral  vanishes. 
If  the  region  of  integration  extends  to  infinity  the  surface  integral  vanishes 
since  the  elastic  displacements  vanish  (at  least  as  fast  as  1/R)  as  | r— r | 
and  |r-r2|  go  to  infinity.  If  the  region  is  finite  we  assume  the  usual  homo- 
geneous boundary  conditions  such  that  either  the  displacement,  u (r) , or  the 

m 

normal  stresses,  T = T n = c .u  . (r)n  , vanish  on  the  free  surface. x 
m ms  s mpsj  p,j  s 

In  either  case,  one  then  obtains 


V.  X U,  x 

uy  — 2 — 1.)  = uv(— 1— 2 


(17) 


Note  that  whereas  the  reciprocity  of  Eq.  (8)  applies  only  to  an  infinite  iso- 
tropic homogeneous  medium,  the  reciprocity  of  Eq.  (17)  is  completely  general. 

B.  The  Integral  Formulation  of  Scattering 

The  simple  solution  of  Eq.  (14)  applies  to  an  infinite  homogeneous, 
isotropic  medium.  Although  such  a solution  has  only  a limited  idealized  appli- 
cation per  se,  it  can  be  used  as  the  basis  for  an  important  extension  of  the 
analytic  description  of  wave  propagation  in  inhomogeneous  elastic  media.  We 
assume  that  the  medium  of  interest  is  isotropic  and  homogeneous  for  R “*■  °°, 
but  that  in  the  vicinity  of  the  origin  there  is  an  inhomogeneity  which  is  the 
cause  of  scattering  of  elastic  waves.  We  denote  the  "unperturbed"  region  by 


The  more  general  condition,  that  everywhere  on  the  boundary  we  have 
au_  ^msns  = which  01  and  S are  arbitrary  (but  fixed)  scalar  parameters 

which  may  depend  on  position  or  on  the  component  m,  also  leads  to  the  vanishing 
of  the  surface  integral  point  by  point.  The  condition  may  even  be  further 
generalized  to  the  extent  that  Ct  or  3 may  be  taken  to  be  symmetric  matrices, 
although  such  generalization  is  probably  physically  artificial,  except  for  the 
the  case  where  a and  8 are  diagonal. 
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p(r)  = P°  + Ap(r) 

c . (r)  = c°  . + Ac 
mspj  — mspj  mspj 

The  equation  for  the  displacement 


(18.1) 

(18.2) 

field  in  the  absence  of  applied  forces 


is 


(c  .u  .) 

' mspj  p,jV 


2 / \ 
poo  u (r)  = 
m — 


This  is  then  written  in  the  form 


(19) 


c°  .u  . (r)  + p°oo^u  (r)  + ( Ac  .u  . ) + Apw^u  (r)  = 0 . (20) 

mspj  p,js  - m — \ mspj  p,j  ',s  ra  — 

When  Eq.  (20)  is  compared  with  Eq.  (3),  it  is  clear  that  the  last  two  terms 

play  the  role  of  an  impressed  force  on  the  medium  (albeit  one  whose  magnitude 

at  each  point  is  dependent  upon  the  displacement  field  at  that  point).  If 

c°  . is  given  by  Eq.  (2)  we  can  wrice  a formal  solution  in  terms  of  the 

Green's  function  for  the  homogeneous  isotropic  medium.  We  are  interested  in 

describing  the  waves  scattered  by  the  inhomogeneity;  we  therefore  express 

u (r)  as  the  sum  of  the  incident  field  u°(r),  which  is  the  wave  field  in  the 
m — m - 

infinite  homogeneous  medium,  and  which  satisfies  the  equation 


o o . . , o 2 o , . 
c . (r)  + p 0J  u (r)  = 0 

mspj  p,js  — m — 


(21) 
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and  a scattered  wave,  (rO  generated  by  the  inhomogeneity, 


u (r)  = u°(r)  + u^  (r) 
m — m — m — 


(22.1) 


umS)<i)  " f ^(£;£>[(4cvspJ(£)“plJ(£>)iS  + ?.C(r,)K2l.,jfr')jjr' 

* A,v  (i-i>  [( ' ‘'vsp 3 & “p , J <I>), , 


+ u^Ap  (r*)u y (r ) dr* 


(22.2) 


The  integral  in  Eq.  (22.2)  in  principle  extends  over  all  space  but  since 

Ac  (jr)  and  Ap(jO  are  such  that  they  vanish  as  | r'|  00  the  integral  in 

vspj 

actuality  is  limited  to  a finite  region  of  space.  The  scattering  cross  sec- 
tion, by  definition,  is  evaluated  in  the  far  field  (i.e.,  at  a position  which 
is  far  removed  from  the  inhomogeneity  which  is  scattering  the  elastic  waves). 
In  defining  the  scattering  cross  section  one  therefore  needs  the  asymptotic 
behavior  of  the  Bessel  functions.  These  are 


, -iz 

h (x)  = 
o z 


(23.1) 


. -iz  r 

h7(x)  = e’: — 3 + 3iz  - z‘ 


(23.2) 


and  hence,  from  Eq.  (14) 


G.  (R)  = 


V fe.  - <2 


. „ X.X  - _ I 

-ikR  , i 2 -ikR  j 

e + — j — k e ] 

R J 
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and  therefore  with,  r = r£2, 


u(s) (r)  = u(s)(rft)  = 7 [”k2(6  - )e~iKrU  (<£) 

m — m — r [_  jm  3 m 3 


+ k2fl.fi  e"lkrU.  (kQ) 
J m J “■ 


+ O 


(?) 


(25) 


or 


u(3)(r)  -i 


Lk2C-0(kE)ae":Lkr  + K2SMJ(ka)*£e‘1<r  + O 


(25.1) 


and 


U.(qft)  = — y /eicj— — •((Ac.  (r)u  (r)) 
3 4TTPO)2  J JSPn  ~ P’n 


4*  0) 


2Ap(r)uj  Or) | dr' 


where  q may  be  either  < or  k.  Then  set  £ = qft  and  write 


Uj(£)  = ^^/pAp(r')u.  (r)  - iqsAcsjpn(r)up>n(r')|  e^^lr'  (26) 


The  strain  field  u.  can  most  easily  be  written  in  terms  of  the  dilatation, 
J » s 

£ = u , and  the  "deviatrix,"  E.  = y (u.  + u . ) - y z6  , which 

s,s  j s 2 j ,s  s , j 3 js 

represents  a pure  shear;  then 
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c.  u = c . E + -r  c . £ 

sjpn  p,n  sjpn  pn  3 sjpp 


(27.1) 


which,  for  an  isotropic  medium,  becomes 


u = 2yE  + (\  + ~ y)e5  . 
i P.n  sj  \ 3 / sj 


(27.2) 


and  X + — y = K is  the  bulk  modulus. 

Equations  (25)  and  (26)  therefore  give  the  scattering  from  a defect  or 
inhomogeneity  in  an  infinite  homogeneous  isotropic  medium  in  terms  of  the 
displacement  and  strain  fields  within  the  inhomogeneity.  To  find  u.  Or ) and 
Up  ^ (r) , however,  it  is  necessary  to  obtain  the  solution  of  the  "internal" 
field  defined  by  Eq.  (22).  Equation  (22.2)  is  first  integrated  by  parts  to 


u (r)  = u(o)  (r)  + u)2  / G . (r-r')Ap(r)u.  (r)dr' 
m — m J m j — — — j ~ ■ — 


/ G . (r-r)  i Ac. 

J m3  ,s ) 3 


(r)E  (r)  + \ Ac  (r')e(r')ldr'  (28) 

spn  — pn  — 3 jspp  — — I — 


C.  Scattering  Reciprocity 

Equation  (17)  defines  a reciprocity  relationship  between  the  displacement 
field  at  point  produced  by  an  impressed  force  at  _r 0 and  the  displacement 
field  at  produced  by  an  impressed  force  ar  r_^.  It  is  also  convenient  to 
have  a reciprocity  condition  which  applies  directly  to  the  scattering  cross 
section  and  the  scattering  vector  defined  in  Eq.  (25). 
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We  consider  a scattering  region  near  the  origin  of  a coordinate  system 
and  two  points  _r^  = r^f2^  and  = *2—2’  assume  that  r^  and  r^  are  both 

large  compared  to  the  dimensions  of  the  scattering  region  and  with  respect 
to  the  wavelengths  of  shear  and  compressional  waves  in  the  medium  (Fig.  1). 

We  describe  the  reciprocity  statement  of  Eq.  (17)  to  be  composed  of  two 
components — the  direct  propagation  of  the  wave  field  from  point  1 to  point  2 
and  a wave  field  scattered  to  point  2 by  the  scattering  region.  Then  we  can 
write 


V . \ 

uy  — 2’— 1 = 


Gy(  l~2_“l I ) + uyS)(^2) 


(29.1) 


Cs ) 

where  is  the  displacement  field  at  which  exists  because  of  the 

presence  of  the  scatterer.  Without  any  loss  of  generality,  one  can  describe 

(s) 

u^  (r^,)  as  the  result  of  a scattering  of  an  incident  wave  by  the  scatterer, 


uyS>  (— 2}  = Ga^rl)Wy(-2  ’~-(1) ,rl) 


(29.2) 


where  we  have  introduced  the  minus  sign  to  indicate  that  the  incident  wave 
is  moving  in  the  direction.  Similarly 


- <(1%-^!)  + G>X(ri;-c(2),r2) 


(30) 


In  writing  Eqs.  (29)  and  (30)  , we  have  done  no  more  than  define  new  quantities 
Wy  (l-2 ’ ri)  aru*  Wv(li!~^(2)  «r2)-  We  are  interested,  however,  in  the 
limiting  case  r^  -*■  °°,  r^  “*■  and  in  this  limit  W^r^ will  be  inde- 
pendent of  the  magnitude  of  r^  and  similarly  W^/r^;-ft^ ,r^ ^ will  be  independent 


15 


Rockwell  International 


Science  Center 


SC579.4FR 


of  r^-  Since  Eq.  (17)  is  valid  for  an  arbitrary  inhomogeneous  elastic  medium, 
it  is  also  valid  for  the  specialization  to  a localized  inhomogeneity.  Further- 
more, the  Green's  function  for  the  infinite  homogeneous,  isotropic  medium  is 
symmetric  so  that  we  have 


Go<ri,"u(V-5ll>-r2)  ■ GS<r2>Uv(v2<2).r2) 


In  the  limit  as  r^  and  r2  become  large,  Eq.  (31)  leads  to 


r„  K2  (5  - 

2 \ va  vc 


= rx  <2(. 


(2)n(2)\  r2  . ,2n(2)n(2)o  2 ,s 


6 6 - e * + k^'OT'e  ‘ tf 


y s y s 


y s 


j(— i*’-— (2)  ,r2) 


Now  W^r^*"^^)  and  VT^_r^ ,r?^  can  be  decomposed  into  longitudinal  and 
transverse  components.  Ar  large  distances  ^ , r ^ j must  behave  as 

_ <(r2:-2(1)-ri) 


r r 

l,r2 


k2fi(2)fi(2)Ua  ( Q(2)  ,£;-S)(1)), 
y s s \ — ’ ’ — / 
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(where  we  have  introduced  the  additional  descriptor  l or  t to  emphasize  the 


longitudinal  or  transverse  components  of  the  vector  U ) , with  a similar 


0 / 2 \ 

expression,  of  course,  for  it  s-fl'  ,r2  • Equation  (32)  can  then  be 


decomposed  into  the  form 


i<(r1+r2) 


k2(6  -n<2>n<2V(o(2),t;-na),t)e  lri 

\ va  v a I [ \ ys  y s ' s\—  — ' 


+ A(2)a(2)U°(!!(2),t;-!!<1),t)e'torikr2} 

y s s \ — / 


-ikr, -i<r„ 


+ kV^Vffi  - ^2>^2>)ua(f2(2>,t;^(1)^)e"  1 

v a I V ys  y s'  s\—  — / 


+ kV2V2>^(2<2>,«:-n(1>i)e"lk<ri+r2)} 


k2(<5  5 - n<2)n<2)){K2(6  -2<1)n<1>)us(n<1>,t;-!)(2),t)a'1'",:i+t2> 

\ y s y s / I \ va  v a I a\-  ’ / 


+ k2a(1>s!aV  n(1),t;-n<2),t)f'lKr2”lkri 

v a \ a — ’ — ’ / 


- n(1)n<1>)D*(n<1),t;-!!<2),j)e'1"1"il'r2 


VC  V a 


+ kV1)^1)U®(n(1),Z;-n(2),Jl)e  ( 1 2>  . (34) 


From  this  we  obtain  the  fundamental  scattering  reciprocity  relations 


1 
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(35.1) 


Ug(fi(2)  jl)  = irj(n(1)  ,s,;-n(2)  ,n) 


(35.2) 


U^(2\t;^1\£)=U^1>5.;-^2>,t) 


(35.3) 


or  simply, 


Ug(il(2);-n(1))  = ;-n(2)j 


(The  fourth  expression  which  can  be  obtained  from  Eq.  (34)  is  equivalent  to 
Eq.  (35.3)  with  s and  a and  (1)  and  (2)  interchanged.) 

The  scattering  amplitudes  defined  in  Eqs.  (29)  and  (30)  and  their 

O' 

decomposition  into  longitudinal  and  transverse  components  defined  in  Eq. 
(33)  refer  specifically  to  the  scattering  associated  with  an  incident  field 
of  unit  amplitude.  The  incident  wave  will  then  carry  an  energy  flux  per 
unit  area  given  by 


nin  3,  in 

Q = pw  /q 


(36.1) 


where  q*n  is  the  propagation  number  (either  k or  <)  of  the  incident  wave. 
The  scattered  energy  flux  per  unit  area  normal  to  the  direction  n is  given 
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= <jdlm 


(a+—  y)eu  + 2yE  .u. 
L'  3 / n nj  j 


(36.2) 


£ * ^ VVl®e'lkr  + ° (7) 


(37.1) 


2E  . = 
nj 


— K3(fi  <5._  + ft. 6 

r L n js  J 


__  - 2fi  S7.fi  )U  (<fi)e 
ns  n j s s — 


-iKr 


2k3(fi  fi  - i 5 ,)fi  U (kfi)e‘lkrl  + O / -U 
' n j 3 nj'  s s - J \r2/ 


(37.2) 


The  scattering  cross  section  is  defined  as  the  far  field  scattered  energy 
flux  per  unit  solid  angle  divided  by  the  incident  plane  wave  energy  flux  per 
unit  area 


S(fi)  = — p-  lim  r2fi  Q 
~ Qin  ^ nn 


in 


lim  r Im 


puj  r-*» 


(x  + \ u)e*fi  u + 2yfi  E* . u . 

\ 3 / mm  mmjj. 


= q 


i 


k3|fi  U (kfi)  I 2 + <3(6  - fi  fi  )ReU*(icfi)U  (<fi)  I 

mm—1  ms  ms  m — s — J 


= qln|k3[fi  U (kfi) | 2 + K3|U(<fi)|2  - <3|fi  U (Kfi)!2 
L mm  — — mm  — . 


(38) 
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S(fi)  = q 


in 


k3 I U (kfi) | 2 + K3|gxU(Kfi) 


SC579.4FR 


(38.1) 


It  is  clear  from  Eq.  (38)  that  the  vector  U(^)  contains  all  of  the 
information  in  the  displacement  field  u at  large  distances  from  the  scattering 
center  and  contains  it  in  a form  which  separates  the  purely  geometrical  aspects 
of  the  scattering  from  the  properties  of  the  scatterer. 

It  is  useful  to  further  decompose  the  scattering  cross  section  into  the 
shear  and  compressional  wave  scattering.  We  then  write 


c /■.■.SC  f\in.  i ^ I /i  i . — . HI . i 2 

& > = k uma(k^  ;_k^  } 


I 


(39.1) 


1 


Vt(-SC;-ln)  = kK3|(6ms  " ^sC)fi0nUsa(<-SC;_k-in)|2 


(39.2) 


W^Vn)  = <k3!C(6va  - ^an>Umv^SC--^n^2 


(39.3) 


St^t(nSC;f)ln) 


k4  I (5  -nscnsc)(6  - (2in.Qin)u  (kflsc;-icnln)|2 

ms  m s va  v a sv  — ’ — 1 


(39.4) 


and  we  obtain  the  reciprocity  theorems  for  scattering 
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(40.1) 

(40.2) 

(40.3) 


D.  The  Bom  Approximation 

The  scattering  cross  section  may  be  calculated  from  Eqs.  (22),  (26)  and 

(38.1).  The  Born  approximation  in  its  simplest  form  assumes  that  the  scattered 

wave,  u (r)  may  be  neglected  in  (26)  and  that  u (r)  may  be  approximated  bv 
m — m — 

the  incident  wave  u°(r).  We  now  write  the  incident  plane  wave  as 
m — 


u"(r)  = b e 
m — m 


. xn 
-i£  -r 


(41.1) 


with  the  understanding  that  one  characterizes  the  longitudinal  and  transverse 
waves  by  the  relationships 


longitudinal : 

in  in  ,2 

q q = k 
nm  nm 

, m , 

, b q = k 

m m 

, b b = 1 

m m 

transverse : 

in  in  2 

q_  q_  = « 

, in  _ 

» b q_  = 0 

, b b = 1 

■m  bn  m m m m 


(41.2) 


(41.3) 


The  scattering  amplitudes  in  the  Born  approximation  are  therefore  given  by 
the  following  expressions  [Cohen,  1976;  Gubernatis,  et  al. , 1976],  when  the 


scattering  inhomogeneity  is  restricted  to  be  isotropic: 
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(a)  Incident  Longitudinal  Wave 


V£  = f COS0  " 


fd  + 2fp(cos20  - |). 


f f 

U„  = — sin0  sin28 

Jl-*-t  p k y 


in  sc 

where  9 is  the  angle  of  scattering,  cosG  = , 


„ „„  1 /.  / \ iK-r.  ir  sc  in 

fp®  = 4^7AP(l)e d£  5 K = £ “£ 


f (K)  = ^ f [AX  (r)  + | Ay(r)[  e1- ^dr 


fy<&  = ^ f^(l)eX--dL 


(42.1) 


(42.2) 


(43.1) 

(43.2) 


(43.3) 


Hence,  these  functions  are  the  Fourier  transforms  of  the  perturbation 
of  the  density  and  elastic  constants  in  the  scatterer. 

(b)  Incident  Transverse  Wave  (Plane  Polarized) 


Jt+£  = cos^ 


f f 

— sin9  sin29 

P < U 


U = cos4> 

t->t,, 


— cos9  - — cos29 

P P 


U = -sin<{> 

t^l 


f 

J 

P 


f 

_V 

y 


— ^ - — cos9 


(44.1) 

(44.2) 

(44.3) 
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where  <J>  is  the  angle  between  the  plane  of  polarization  and  the 
scattering  plane,  and  U , U are  the  scattering  components 

t t(|  t tj_ 

parallel  and  perpendicular,  respectively,  to  the  scattering  plane, 
(c)  Incident  Transverse  Wave  (Circular  Polarization) 


U = I e-ir<^ 

t -+1  2 

r 


f f 

— sin0  — sin20 

.P  K y 


(45.1) 


U = i e"ir(f 

t -*-t  2 

r s 


(cos©  + rs)  — - - (cos20  + rscosb)  — 1 — 

p y 


(45.2) 


1 — ircj) 

= — O 


e (cos0  + rs) 


f f ■ 

— - (2cos0  - rs)  — 

p y 


where  r = -1  gives  the  polarization  of  the  incident  wave  and  s = ±1 
that  of  the  scattered  wave.  We  note  that  for  forward  scattering 
there  is  no  depolarization  (s  = r,  only)  while  for  backward  scattering 
there  is  total  depolarization  (s  = -r,  only). 

To  extract  the  scattering  cross  section  from  the  scattering  amplitude 
one  has 

S («)  = q±nqSC | qSCU| 2 (46) 

4- 

E.  Improved  Born  Approximation 

The  expressions  given  in  the  previous  section  are  valid  only  for  a 
limited  range  of  parameters.  In  general,  one  must  restrict  their  use  to 
cases  not  only  for  which  Ka  is  small  (where  a is  a characteristic  dimension 
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of  the  inhomogeneity)  but  also  for  which  the  magnitude  of  the  perturbation, 

|Ap(r)/|p,  |Ac(_r)/c|  is  also  small.  This  latter  constraint  is  the  more 

serious  one  since  it  means  that  the  Born  approximation  is  inapplicable  to 

most  problems  of  interest  to  NDE,  where  the  most  common  scattering  center 

is  a void,  for  which  Ap  = -p.  The  Born  approximation,  in  replacing  u (r) 

m — 

by  u°(r)  assumes  that  the  integrals  in  Eq.  (28)  may  be  neglected.  This  is 
m — 

an  acceptable  approximation  for  u (r)  but  it  is  not  acceptable  for  its 

m 

derivative  or  for  £ . (r) . To  see  this  we  note  that  u (r)  will  involve 
mj  — ra,s  — 

G . (r)  which,  as  may  be  seen  directly  from  Eq . (5),  must  contain  delta- 

mj  ,ps  - 

functions . 

Specifically,  we  find,  as  extensions  of  Eq.  (14) 


VR) 


[A  (kR)6  . - 6.  4>.(o)(R)  + x , x <t)^o)  (R)"j 

ATTipco2  L 0 1 Jm  2 - 


where 


4,<P)(R>  - i Kn+p+1h  <*R>  - k"+I,+1h  (kR) 
n _n  [_  n n 

K 


From  Eqs.  (12)  and  (13)  one  finds 


*(p)(R)  = rVp}(r)  + o(p+2)(r) 

n 2n+l  L n+1  n-1 


(49.1) 


- \ j"(2n-l)<}>(p)  (R)  - 4>(p^2)(R) 
R2  L n-2  J 


(49.2) 


3$(P)(R)  /_n  (p) 

_5_ = ^P>(R)  = -xo<D  P 
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where  the  delta  function  is  defined  such  that 


R 6(R)  = 0 


nf( 


4 it  / 6 (R)R  dR  = 1 


The  A and  X are  totally  symmetric  expressions,  [Cohen,  1976] 
j ms  p j ms  p 


A.  =5.6  +5.6  +6.6  (52.1) 

jmsp  jm  sp  js  mp  jp  sm 


X.  = x.x  6 + x x 6.  + x.x  6 + x x 5.  + x.x  5 + x x 6. 

jmsp  j m sp  spjm  j s mp  mpjs  jpms  msjp 

(52.2) 


with  the  properties 


A.  = A.  = 56.  , A 15 

jm  jmss  jm  ss 


X.  = X.  = R25.  + 7x.x 

jm  jmss  jm  j m 


X = 10R2 
ss 


(52.3) 

(52.4) 


Therefore,  the  first  curly  braces  in  Eq.  (51.2)  vanishes  if  contracted. 
Since  it  may  be  verified  that,  for  R 0, 


*(P)(R)  , <2  - ^ 
n (R)  „2n-l 


p = 0 


<P  - kP 
„2n+l 


p * 0 


-2  - i 

it  is  seen  that  G.  (R)  R and  G ~ R 

jni,s  jm,sp 
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Equation  (28)  may  be  differentiated  to  obtain 


Vp®  ■ u»!p<*>  +“2/G„j, 

+ /Gmj , sp  i * A (I‘J ''  j s £ (-r-‘  + 24p(r')ej8(r')fdr'  . 


(52) 


We  define  symmetrized  functions 


G,  (R)  = G . . (R)  = ir 
J >®P  mp;j  2 


G . (R)  + G . (R) 

L ®J  ,p  pj  ,m 


2 

4TTipw 


" To  (KR) (x  5 . + x 5 . ) 

2R  1 p mj  m pj 


+ 


(x . 6 + x 6 . + 

j mp  m jp 


x 6 . )<j>^o)  (R) 

p jm  2 


. (o) 

x.x  X $ 
j m p 3 


(53) 

and,  with  the  delta-function  component  omitted, 


G 

nip;  js 


, G . (R)  + G . (R) 

L mp;j ,s  rap ;s ,j 


1 

4iripu)' 


h (kR) (6  5 . 

2R  1 ps  jm 


+ 5 .5  ) 

PJ  ms 


+ 7 h (kR)(o.  xx  +6. xx  +6  x.x  +d  xx) 

1 jmps  jpms  ms  j p psmj 


- * 


4>.J0)(R) 
jmsp  3 


+ x.x  x x (R) ; 
j m s p 4 


(54) 
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We  now  introduce  an  operator  notation  which  can  be  defined  in  an  obvious 
way  by 


£ ,u  = oj2  / G . (r-r‘)Ap(r‘)u.(r,)dr' 
®J  j J mj — J — 


(55.1) 


fC  c 

■ m;  js'js 


- f G . (r-r*)Ac  . (r*)e.  (r^dr* 
J m;pn pnjs  — js  — — 


(55.2) 


mp 


5 jUj  " 


(55. ”) 


'Jmp ; j s£  j s Ap ; nq  Acnq  j s ^ £ j s ^ d- 


(55.4) 


Then  Eqs.  (28)  and  (52)  lead  to 


u (r)  = u“(r)  + 
ra  — m — 


<g  .u. 
®J  J 


rQ  . e. 

m;  js  js 


(56.1) 


e (r)  = e"  (r)  + 

mp  — mp  — 


fg  . u . + fQ  . £ . 
mp;j  J mp? js  js 


15poo 


(<2  - k2)6  Ac 

mp  nn j s 


(3k2  + 2k2)Ac 


■ e,  (r) 
mpjsJ  js  - 

(56.2) 
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where  the  final  term  in  Eq.  (56.2)  arises  from  the  delta  function  in  Eq.  (51). 

Introducing  the  shorthand  notation  'Q  = rS  ,u.  + . e.  for  the 

mp  mp ; j j mpjjs  js 

integral  terms  in  Eq.  (56.2)  gives,  formally,  an  algebriac  matrix  equation  to 

be  solved  for  £ (r) . However,  in  the  case  where  Ac  . is  isotropic  this 

mp  — mpjs 

system  can  be  diagonalized.  In  the  isotropic  case,  which  is  the  only  one  we 
shall  explore  in  detail  from  this  point  on,  Eq.  (56.2)  is 


0 ' ' • L 


£ n (r ) = £ (r)  + a 

mp  — mp 


mp  ISpc2 


5<2AX(r)e(r)6  - 2 (<2  - k2)Ay  (r)£  (r)6 

~ ~ mp  — — mp 


+ 2 (3<2  + 2k2)Ay(r)£  (r) 

— mp  — . 


(57) 


Evaluating  the  dilatation  and  deviatrix  of  this  expression  gives 


. AX(r)  + ~ Ay(r) 

£ (r)  = £ (r)  + (Sss  ~ a(r)£(r)  ; a(r)  = ZT~2u 


(58.1) 


E (r)  = E°  (r)  + 'g  - | " b(r)E  (r) 

mp  — mp  — mp  3 ss  mp  mp  — 


, , x 2(3A+8u)  a , \ 
b(r)  = 15yTA+2y)  Ap(r) 


(58.2) 


and  hence 


f.  (r)  = 

mp  — 1+b  (r) 


L<p<r)  + ®.P  + +«..) 


a(r) 


b (r)  - a (r) 
1 + a(r) 


(59) 
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The  integral  equation  (59)  is  not  yet  in  a convenient  form  for  a Neumann 

expansion  since  'Q  still  involves  integrals  which  do  not  vanish  in  the  limit 
mp 

as  the  characteristic  dimension  of  the  scatter  goes  to  zero,  although  for  the 
special  case  of  a sphere  (because  of  its  symmetry)  those  integrals  vanish. 

In  that  one  case,  then,  using  £ = (1+b)  1 (c°  (r)  + ^ ae°(r)6  \ in  Eq.  (26) 

mp  \ mp  — J — mp  / 

gives  the  "corrected"  Born  approximation  [Cohen,  1976;  Kal  and  Knopoff,  1976]. 

In  order  to  evaluate  the  complete  first  order  approximation  to  the  strain 

field  inside  the  scatterer  it  is  necessary  to  evaluate  the  lowest  order 

contributions  of  the  integrals  fSm^  > and  hence,  equivalently,  the  long-wavelength 

limit  (static)  Green's  function.  If  we  assume  the  long-wavelength  limit,  then 

u (r)  will  be  a slowly  varying  function  of  position  within  the  scatterer  and 

we  may  use  a mean-value  theorem  to  remove  u (r)  from  the  integrals.  Further- 

m — 

more,  we  may  use  Eshelby's  result  that  e is  constant  inside  a homogeneous 
ellipsoid  in  the  static  limit  and  remove  it  from  the  integral  in  that  case 
without  resort  to  any  approximation.  Then,  keeping  only  those  terms  which 
can  contribute  in  lowest  order,  we  may  write  Eqs.  (58.1)  and  (58.2)  in  a con- 
venient symmetric  form: 

(l+a)e(r)  = e°(r)  + (aA  + \ Au)ge(r)  + 2Ayg  E.  (r)  (60.1) 

(l+b)Emp (r)  - E°p(r>  + (aA  + f Au)gmpE(r>  + 2Aygmp ;J sEjs<r)  (60.2) 

where  the  integrals  are  defined  by 
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= / C ..(r-r)dr 
J pp; j J — - 


(61.1) 


g = f G . . (r-r')dr'  - i 

mp  J mp;jj — 3 


g<5 


mp 


(61.2) 


mp 


= / G . (r-r^dr’ 
;js  J mp;  js — 


(61.3) 


We  note  that  it  is  possible  to  write  g.  in  Eq.  (60.1)  as  the  same  functional 

J s 

expression  as  g in  Eq.  (60.2)  because  of  the  fact  that  E,  has  zero  trace, 
mp  j s 

Therefore,  adding  the  term  - -r  5.  G (r-r)E.  (r)  in  Eq.  (58.1)  adds  nothing, 

3jspp;nn js  — 

but  allows  the  symmetric  structure  of  the  equations.  In  fact,  since 
_ 22  — ikR 

G (R)  = -(k"/47rp<jj  R)e  is  a lower  order  singularity,  its  contribution 

pp ; nn 

2 2 

to  g.  is  (Xk  a ) and  hence  vanishes  to  second  order.  One  must  be  cautious 
1 s 

here;  G . (R)  still  contains  a delta  function,  -6 (R) / (X+2u) , which  is,  of 

PP  > nn 

course,  not  a lower  order  singularity  and  cannot  be  neglected. 

The  integrals  must  now  be  evaluated  in  the  long-wavelength  limit. 

Expanding  the  Bessel  function  gives  us 


jm 


4iTpw2R 


y (2k2  + k2)<5 , - 5.  W 

3 !m  \ p2  jm/ 


2 

8ttpcj  R 


4,-^) 


(62.1) 
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or,  in  the  long-wavelength  limit. 


jm  8TTy(A+2y)R 


x.x 


(A+3y)6. ■ + (A+y)  -Li" 

jm 


R J 


(62.2) 


G.  (R)  = =-=• 

ism  „ 2„3 

J 87rpco  R 


2 

< x 


X.X  X 

3 - - x 6 -x  5 ,-x  6. 

„2  j ms  s mj  in  js 


>- 


(63.1) 


G° 

js;m 


-1 


8fry  (A+2y)R“ 


x.x 


y(6  x +6  x.)  + (A+y)x  3 — 1 V5- 

jm  s sin  j vi-/  mi  ^2 


5. 

J s . 


(63.2) 


and 


Gmp;js(R)  = 23 

l6irpco  R 


„ ( ~ 30x  x x.x 

ic  .26  6.  - — «■  (x  x 6 . + x.x  6 + X . ) + ) — — 

| mp  js  r2  m p js  j s mp  mpjs  R4 


o | X 15x  x x . x | 

-2k2  A . -3^EJs+  + 

I mPJs  r2  r4  I 


(64.1) 


G°  . 
mp;  js 


-1 


167ty  (A+2y)R“ 


o..a  , 3A  „ 30(A+y) 

X . - x x x.x 

rapjs  r2  mpjs  r4  mpjs 


)+(,3^.4  v •' 

I mP  \ R2  is/  \ r2  mp/  jsl_ 


(64.2) 
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In  these  expressions  only  the  lowest  order  terms  have  been  retained;  the 
Green's  function  is,  to  this  order,  the  Green's  function  for  the  elastostatic 
problem  [Eshelby,  1961;  Love,  1934]. 

We  now  use  the  divergence  theorem  to  convert  these  integrals  into  surface 
integrals  and  evaluate  them  at  the  center  of  the  ellipsoid.  This  transforma- 
tion, however,  recaptures  the  effect  of  the  delta  function  in  G . and  it 

mp ; j s 

will  therefore  be  necessary  to  consider  a small  sphere  which  excludes  the 
origin  as  part  of  the  surface  of  integration  or  *0  appropriately  modify  the 
equation  in  order  to  avoid  including  this  contribution  twice.  We  shall  choose 
the  latter  approach  and  merely  insert  the  modifications  later.  Then 


x^  = a^sin0cos<}) 

2 2/2  2 2 2 0 o 

= a^sinOsin^  , R = a^  + ( a^  - a^  + (a0  - a^)  sin“d))  sin“9  (66) 

x^  = a^cosS 
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ala2a3 

n.dS  = 2 — x sin0d8d<}>  , 

J 


(67) 


and  g . can  be  expressed  in  terms  of  two  classes  of  integrals 
mp ; js 


1/3 


mp 


4tt 


/ 


x x sinOdSdd) 
m p 


(68.1) 


\l/3 


K 


mpjs 


^a^a^a^)  ' J J"  x_x_xJx^sin0d6dcf> 


m p j s 


4tt 


(68.2) 


From  symmetry  we  find  immediately  that  J is  diagonal  and  that  there  are 

mp 

only  six  independent  non-zero  components  for  K jg.  Using  six-vector  notation 

with  mp  a,  is  ^ B we  see  that  K . may  be  written 
r J mpjs  7 


K . = K 0 

mpj  s aS 


K1  K6  S 

K6  K2  K4 

k5  k4  k3 


K.  0 0 

4 

0 K5  0 

0 0 K, 


(69) 


and  therefore 
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8mp ; j s 


I \2/3 

(ala2a3) 


4y (A+2y) 


(X+y)(6  J. 

mp  js 


- 3K  . ) 


mpjs 


J-+X 
2 2 
a . a 
3 s' 


PL. 


+ 6 


ms 


J . 


6 -JS3-  + 6 -BJ- 


jm  a a pj  a a ps  a a ms  a a.  ] 

Jps  J m s m j PJ/J 


where 


g, 


.( 


\ 2 / 3 

ai  2a3/  r 


otS  2y(A+2y) 


+ UBafJ 


3K2-J22 


3VJ33 


3K4-J22 


3VJ33 


3K4(a2+a3) 


2a2a2 

2a3 


3K5(a^-a*) 


, 2 2 

2ala2 


(70.1) 
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Dag 


0 

J22/a2 


0 

0 

J33/a3 


Jlla3+J33a2 

/ 2 2 
4a2a3 


Jlla3+J33ai 

4a2a2 

4ala3 


Jlia2+J22ai 


, 2 2 
4aia2 


/ 


(70.3) 


Among  the  nine  integrals  there  are  three  immediate  relationships 


J 


J 


J 


11 

22 

33 


K,  + K,  + Kc 
lop 

K6  + K2  + K4 
k5  + k4  + k3 


(71.1) 


Since  it  is  also  clear  from  Eas.  (61.3)  and  (64.2)  that  A . is  a symmetric 

ab 

matrix,  the  asymmetric  form  of  Eq.  (70.2)  must  be  illusory  and  it  can  indeed  be 
verified  that 


3 (a2  - 

a3)iC4 

' a2J33 

- a2J 
a3  22 

3 (a3  - 

al)K5 

' a3JU 

- ^33 

(71.2) 

2 

2, 

2 

2 

3(ax  - 

VK6 

= alJ22 

" a2Jll 
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so  Chat  the  six  may  all  be  expressed  in  terms  of  the  three  J's. 
One  also  finds 


8mp  8mp ; j j 


i+2y)  /(xmn 


87r(A+2y) 


% ds 

+ x n ) — =- 
P P m r3 


1 (ala2a3)‘ 


-i-  \ x j • j — - <5 

X+2y  a a mp  3 mp 

m p 


where  the  second  term  in  the  bracket  comes  from  subtracting  out  the  delta 

function  in  G ...  We  may  then  write  Eq.  (60)  as* 
mp;jj  H 


E = E / (1  - 4Ayg  ) 

a a 6aa 


a = 4,5,6 


(73.1) 


i?  + 2 Ay  ! i 


1 /Va 


E!  E1  + 2Ay  |(8H  ~ 813)E1  + (g12  " 813)E2  “ \+2y  ^ 2 j J11  ~ 3 £ 


(73.2) 


E°  + 2Ayj(g12 


823)E1  + (g22 


823)E2 


(73.3) 


It  should  be  remembered  that  in  6-vector  space  a scalar  product  is  defined 
by  A- B = A B + A_B„  + A.B,  + 2(A.B.  + ACBC  + A, B. ) in  order  to  account  for 
the  svmmetrization  of  the  subscripts  ij  and  ji  into  a single  component. 
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/‘■'ll 

(l+a)e  = e°  - AKc  [ —5— 


(73.4) 


2/3  2 

where  c = (a^a^a^)  /(X+2p)  and  AK  = AX  + — Ay,  and  use  has  been  made 

explicitly  of  the  trace  condition 


E±  + E2  + E3  = 0 


For  a spherical  scatterer,  a^  = = a and  = J ^ = -J33  = 

K-,  = K„  = K,  = K = Kc  = K,  = vV.  Then  we  find,  from  Eq.  (70)  that 
1 Z J j 4 j DID 

-2Ay(gi;L  - g13)  = -2Ay(g22  - g23)  is  just  the  "shielding"  factor,  b,  previously 
obtained  from  the  delta  function  in  Eq.  (58.2)  and  that  this  corresponds  also 
to  -4Ayg^  = -4Ayg33  = -4Ayg^^  = b.  Thus,  for  a spherical  scatterer,  the  only 
contributions  to  the  renormalization  of  the  internal  strain  field  are  the 
shielding  factors  a and  b which  arise  from  the  delta  function  singularity  in 
Green's  function.  For  an  arbitrary  ellipsoid,  or  for  any  more  general  shape, 
we  have  now  explicitly  exhibited  the  volume  "polarization"  contributions  to 
this  renormalization. 


Ellipsoid  of  Revolution 


For  an  ellipsoidal  scatterer  with  a^  = a2  we  then  have  J = J 22 , 
Kx  = K2  = 3K&,  K4  = K5,  gu  = g22,  and  g13  = g2y  We  also  find 


w 
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- [a3  - • 


11  22 


/a3  - 0 


-1 

sec  u 

/2  7 

v u “1 


u > 1 (75.1) 


I(u)  = 


, (Va3 )2/34  rTtl  ? . ,] 

“ (.*  - a3)  1 ‘ 1 3 J 


sech  Lu 


l ATS 


u < 1 (75.2) 


K4  = K5  = 


T 2 T 2 

Jlla3  J33al 


(al/a3) 


2/3  2 


H - ai) 


^ r 2 

2 \2~  [a3  + 2al 


2 , 2 


3a1I(a1/a3) 


(75.3) 


K1  = K2  = 3K6  = 4 (J11  - V 


/ \2^3 

/al) 

2 212  L' 

-al) 


iai/a3 ) 
ala 


2a 3 - 5a2a2  + 33^1(3^3^ 


(75.4) 


K3  " J33 


(ai/a3)2/3  r 2 

" 2K4  = 3-p  2\2  L3a3I('al/a3) 


I 2 2V 
(a3  - al) 


,2  2 

4a3  + al- 


(75.5) 


Then,  from  Eq.  (70.2) 


( ai/a3 )2/3  - 

*11  * *22  = '(X-U)JH  - 3(X+y)Kl- 


(76.1) 


|a3/a1)2//3(^+U) 


*12  = 2u(X+2y)  (J11  " V 


(76.2) 
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(76.3) 


(a./a,)2'3  . 

*33 " -Wrlirr  [<X-“)J33  - 3(X+1J)K3] 


(76.4) 


g44  = g55 


2y(X+2y) 


: \ 


yJ„  + 3(X+y)K.  + 


jyJ33  + 3(X+y)K4J 


(76.5) 


(V‘i)2/3  r 

*66  ■ ahx+ffiT  Kl  + 3(X+*‘,K6. 


(76.6) 


Using  these  expressions  reduces  Eq.  (73)  to  a simpler  form,  which  may  be 
solved  rather  directly  to  give 


£°(1  + Y+)  - YcAK(E°  + E°) 
(1+a) (1  + y+)  - 4BycAyAK 


(77.1) 


where 


ill  i33 
2 ~ 2 
al  a3 


(V*!)' 'Xl 


Y = 2Ay(2g13  - glx  - g12) 
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(77.2) 


(77.3) 


These  equations  then  provide  the  shielding  factors  for  the  extended  Born 

approximation.  Given  the  incident  strain  field  £°  (r)  = E°  (r)  +46  e°(r) 

pm  — mp  3 rap  — 

we  use  Eqs.  (73.1)  and  (77 . 1) - (77 . 3)  to  provide  an  accurate  first  approxi- 
mation to  insert  into  Eqs.  (25)  and  (26)  for  the  scattered-wave  amplitude 
or  Eqs.  (26)  and  (38)  for  the  differential  scattering  cross  section,  noting 
however  that  the  components  of  the  displacement  and  the  strain  considered 
here  are  specifically  written  with  respect  to  a coordinate  system 
corresponding  to  the  principal  axes  of  the  ellipsoid. 

The  important  feature  of  this  solution  is  that,  to  this  level  of 
approximation,  the  general  structure  of  the  scattering  vector  U given  in 
Eqs.  (42)- (45)  above  is  unchanged.  The  density  component,  that  involving 
Ap/p  is  exactly  as  given  in  those  equations.  The  dilatation  and  shear 
components,  involving  f^  and  f^  will  still  involve  the  Fourier  transform 
of  the  scattering  volume  and  the  shielding  factors  given  in  Eqs.  (77.1)- 
(77.3).  In  addition,  the  purely  geometrical  aspects  of  the  scattering  law, 
(those  due  to  the  projection  operators  (2  and  6 - (2  Q of  the  incident 

and  scattered  wave  directions)  will  contribute  no  more  than  an  "isotropic" 
component  and  components  proportional  to  the  second  order  spherical 
harmonics,  and  therefore  no  more  than  terms  of  the  form  a^  + a?cos29  or 
sin29  in  the  scattering  angle,  0.  The  amplitudes  of  the  various  components 
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will,  however,  also  depend  on  the  orientation  of  the  scattering  ellipsoid 
with  respect  to  the  plane  of  scattering. 

The  full  implementation  of  this  stage  of  the  analysis  has  not  been 
completed. 


F.  Variational  Theorems 

Although  integral  equations  are  often  much  more  difficult  to  solve  than 
the  corresponding  differential  equations,  they  are  often  the  more  convenient 
starting  point  for  approximate  solutions.  This  is  particularly  true  when  a 
complete  solution  is  not  needed.  In  the  scattering  of  elastic  waves  one  has 
an  ideal  example  of  this  sort  of  approximation.  We  are  not  interested  in  a 
complete  solution  of  Eq.  (3)  or  even  a complete  solution,  for  all  ir,  of  Eq. 
(28).  What  we  are  interested  in  is  the  vector  U(q)  which  is  given  in  Eq.  (26) 
as  an  integral  over  the  volume  of  the  scatterer.  This  situation  lends  itself 
ideally  to  attack  by  variational  methods.  To  establish  a methodology  for  such 
an  approach*  we  consider  an  unknown  function  u(x)  as  the  solution  of  the 
integral  equation 


u(x) 


u°(x) 


/H  (x  ; x’)  u (x)  dx' 


(78) 


and  for  which  we  want  to  evaluate  the  integral 


J 


(x)u)x)dx 


(79) 


In  order  to  define  a variational  approximation  for  J we  introduce  an 
adjoint  equation 


A different  approach  to  the  formulation  of  a variation  equation  has  been  given 
recently  by  G.  S.  Kino  (Kino,  1976]. 


42 


p 


Rockwell  International 

Science  Center 


SC579.4FR 


v (x) 


f(x) 


(80) 


where  H+(x;x)  is  the  adjoint  operator,  such  that  for  any  two  arbitrary 
sufficiently  continuous  functions,  f(x)  and  g(x). 


JJg(x)}l(x;x)£  (x')dxdx'  = If (x'jxJgCxjdxdx^ 


When  H(x;x)  involves  only  algebraic  factors  but  no  operators  (such  as  3/3x, 
3/3xf,  etc.).  H+(x,;x)  = H(x;x;). 

Now  if  we  write 


Q[u,v] 


f‘ 


(x)v(x)dx  + / u(x)f(x)dx 


- j" u(x)v(x)dx  + J'J' v(x)H(x;x,)u(x,)dxdx’ 


(81) 


we  have  a variational  formulation  for  the  integral  J.  If  the  function  u is  equal 
to  u,  the  solution  of  Eq.  (78),  we  find 


Q[u,v] 


/ 


u(x)f(x)dx  = J 


(82.1) 


and  if  v is  equal  to  v,  the  solution  of  Eq.  (80)  , we  find 
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u,v]  = fu 


(x)v(x)dx 


= J'u(x)w(x)dx  - J'J'[{(x;x)\i(x)v(,x)dxdx 


= J f (x)u 


(x)dx  = J 
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(82.2) 


Thus,  if  u and  v are  approximations  to  the  exact  solutions  of  Eqs.  (78)  and 
(80),  Q[u,v]  is  a second  order  accurate  approximation  to  J.  If  we  set  u = u + n» 
v = v + £ we  find 


Q[u+n,v+£]  = J H(x;xf)ri  (x^dxdx* 


-/ 


C(x)n(x)dx 


(83) 


so  that  the  error  is  dependent  on  the  product  of  the  errors  of  the  approximating 
functions. 

The  variational  expression  Eq.  (81)  may  be  easily  improved  for  a scalar 
function.  Let  us  write  u = c^u^,  v = c,,^  where  and  c 2 are  constants.  Thus 
if  we  can  assume  a shape  for  u and  v we  may  use  the  variational  process  to 
define  the  optimum  amplitudes  for  those  functions.  Then,  with  an  obvious  short- 
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We  now  differentiate  Q with  respect  to  and  c2  and  set  these  expressions  to 


(u^f)  = c2(u1,v1)  - c2  (u^Hv^) 


(u  .v^  = ^(u^v^  “ ci^uiHvi^ 


and  from  this  we  find,  [Goldstein  and  Cohen,  1962] 


(u1,f)(u  , ) 

J(U1'V  ■ (u1,vl)  - (v^) 


Equation  (85)  then  represents  a formulation  of  the  variational  problem  which 
is  homogeneous  of  degree  zero  in  the  trial  functions  and  hence  independent  of 
their  normalization. 

To  apply  this  formalism  to  elastic  scattering  we  write  Eqs.  (28)  and  (26) 

in  terms  of  a nine-component  vector,  z = (u  ,£  ,E, ,E„,e)»m  = 1,2,3,  a = 4,5,6, 

— m a i 2. 

with,  of  course,  E^  = -E^  - E2>  and  £g  = -j  £ + Eg,  8 = 1,2,3.  Equations  (26) 
and  (28)  are  then  converted  into  expressions  of  the  form 


uj(£)  = y*fs(r;£,j)zg(r)dr 


(86.1) 


zs^> 


(86.2) 
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z°(r) 


O , . 

£a(x) 


E°(r) 


E°(r) 


o , . 
£ (r) 


(86.3) 


The  expressions  for  the  matrix  elements  H can  be  written  out  explicitly  by 

comparing  the  form  of  (86.2)  with  Eqs.  (55)— (59) . Although  this  identification 

is  straightforward,  the  explicit  expression  for  9X9  matrix  is  cumbersome  and 

need  not  be  given  here.  The  vector  fg(r_;£,j)  may  also  be  found  from  Eq.  (26). 

For  f we  find 
s 


f s D 


e^  — / Ap  . 2iAy  2iAp 

4tt  l o ’°’0’0’  3 q3  ’ 2 q2  ’ 

\ pco  pco 


2iAu  n 

2 ql’°* 

pco 


iAK  \ 

~2  qi/ 

PCO  / 


(87.1) 


fs<>;a,2) 


°,  f ,0, 


q3>°. 


2iAu  n 

T V0’ 

PCO 


(87.2) 


c f . ^ _ e ^ — L n Ap  2lAp 

fs(£’£»3)  - 4Tr  l0»°>  p » 2 

\ pco 


2 q2  * 


2iAy  „ 2iAu  2iAu  iA» 

— qx >o , — q3,  — q3.  - — 

pco  pco  PCO  pco 


(87.3) 


The  generalization  of  Eq.  (85)  is  straightforward.  One  introduces  two 


trial  functions  (vectors) , u (r)  and  v (r)  from  which  one  constructs  two  vectors 

s — s — 

W and  W+  with  components 
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w = f f (r;£,j)u  (r)dr  (88.1) 

s J s — s — — 

(not  summed  on  s or  p) 

w+  = f z°(r)v  (r)dr  (88.2) 

P y P - p — ~ 


and  a 9*9  matrix  with  components 


Wnm  = / V_(Ou  (r)dr5  ” 11  ^ Hnm  ^ dr'dr 

pm  J p — m — — pm  J J p — pm m 


Then  the  variational  expression  for  IF  becomes 


U.[u,v]  = WW-1W+ 
1 o 


(88.3) 


The  variational  solution  for  the  scattering  function  U (^)  is  to  be 
found  by  introducing  an  appropriate  set  of  free  parameters  into  the  trial 
functions  u and  v and  to  determine  the  "best"  values  of  these  free  parameters 
by  requiring  that  IL  be  stationary  with  respect  to  variation  of  those  parameters. 
In  fact,  since  Eq.  (89)  is  already  based  on  an  optimization  procedure  with 
respect  to  a set  of  linear  parameters,  a significant  improvement  in  the  evalua- 
tion of  Uj  comes  from  inserting  an  appropriate  first  approximation  for  u and  v 
into  Eq.  (89)  and  evaluating  the  expression  with  no  additional  free  parameters. 
Thus  we  can  use  ug ( r)  = z°(r)  and  vg(_r)  = fg(r)  as  a first  approximation.  The 
integrals  involved  are  essentially  just  those  required  for  the  Born  approxima- 
tion (i.e.,  Fourier  transforms  with  respect  to  the  change  in  propagation  vector). 
If  the  wavelengths  are  not  too  short  with  respect  to  the  dimensions  of  the 
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scatterer  the  Green's  function  in  H (r-r)  in  Eq.  (89)  may  be  expanded  in 
powers  of  the  frequency  squared.  The  required  integrals  may  then  be  evaluated 
in  a straightforward  manner. 

When  one  uses  u (r)  = z°(r)  and  v (r)  = f (r)  one  obtains  the  w (r)  = w+(r) 
s — s — s — s — s — s — 

which  allows  the  "corrected"  Born  approximation  for  the  differential  scattering 
cross  section,  Eq.  (-38.1),  to  be  written  in  the  form 


S4(«) 


in. 
q k 


3 s=l 


* sc  . . 

(a  >j) 


(90.1) 


and 


st«p 


in  2 
9 


E E ws(£sc.j) 

j s 


E 9jC  E ws(£SC.j) 


(90.2) 


The  lowest  order  variational  calculation  leads,  for  example,  instead  of 
Eq.  (90.1),  to  the  expression 


in 

q k 


7 

Escr>  . sc  . 

9,  2-j  w(q 

j J S=1  S 


)4> 


(91.1) 


where 


9 t 

$ = y (w-i\  w 

s \ o / sp  p 


P=1 


(91.2) 
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In  the  lowest  order,  long-wavelength,  approximation  to  W , the  correction, 

— o 

4>s , is  exactly  equivalent  to  the  shielding  factors  calculated  above  in  Eq.  (77). 

The  accuracy  of  this  variational  formulation  when  applied  to  ellipsoidal 
scatterers  can  be  judged  by  comparing  the  variational  calculation  for  a 
spherical  scatterer  with  the  available  exact  solution. 
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